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1 Introduction 



In recent years, since the discovery of high-T c superconductivity ||, there has been 
a strong increase of interest in low- dimensional strongly correlated electron systems. 
The reason for this is that there are indications that two-dimensional system may share 
certain features with their one- dimensional analogs 0. In one dimension, the Bethe 
Ansatz technique can allow one to exactly solve Hamiltonians in special cases. The 
prototypical examples of such systems are the Hubbard || and t—J models at its super- 
symmetric point 0, |5| as well as their generalizations |], [7], ^ [J . Recently many other 

Among these an 



correlated electron models have been formulated |TjJ, [ITj, [12], [13], [1 
interesting subclass corresponds to isotropic and anisotropic integrable models asso 



dated to the Temperley-Lieb algebra [15, 16, 17], which were shown to be quantum 



group invariant for special choices of the boundary conditions fL8| . Moreover, the 
master equation governing the dynamics of simple diffusion and certain chemical reac- 
tion process in one dimension gives Hamiltonians which are realizations of Hecke and 



Temperley-Lieb algebras |]19| . So, increasing the interest in study quantum Hamiltoni- 
ans derived from representations of the Temperley-Lieb algebra. 

In this paper we will present and solve via coordinate Bethe Ansatz a class of quan- 
tum Hamiltonians derived from representations of the graded Temperley-Lieb algebra. 
The spectra of these Hamiltonians are obtained using various boundary conditions 
(open and closed). We believe that the Bethe Ansatz technique used is general and 
may be utilized in all one-dimensional strongly electronic models associated with the 
graded Temperley-Lieb algebra. 



The Temperley-Lieb (TL) algebra f2(| has been widely used in the construction of 



Yang-Baxter equation solutions J21], [23| , which is a sufficient condition for integra- 
bility of one-dimensional models. It is a unital algebra generated by the operators U k 
, k = 1,2, N — 1 subject to the following constraint 



ul 

U k U k ±iU k 



{Q + Q- l )U k 
U t U k , \k-l\>l 



with Q e C a given number. 
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The connection of TL algebra with the quantum Yang-Baxter equations is well 



known. In |23[ the graded representations of TL algebra were used by Zhang to con- 



struct solutions of the graded Yang-Baxter equation 

(/ ® R(u))(R(u + v) <g> I) (I <g> R(v)) = (R(v) <g> /)(/ ® + v))(R(u) <g> /) (1.2) 

where J G End is the identity operator, and i? G End (V <g> V), with V = Vq © V\ 
a Z 2 -graded vector space. The operator R(u) is obtained from graded representations 
the TL algebra as: 

„ , . sinhfr] — u) smh(u) TT _ „ ,„ „. 

i2 fc (u) = ^ -+ . ^ ! U k , k = l,2,...,N-l 1.3 

sinn r] sinn 77 

where u G C is the spectral parameter and the parameter 77 is chosen so that 

2 cosh r/ = g + g- 1 (1.4) 

The i?-matrix ( |1.3|) gives rise to a solvable vertex model on a planar square lattice. 
The favorite method to apply in such a case is the algebraic Bethe Ansatz. Here one 
uses the fact that the Yang-Baxter equations can be recast in the form of commutation 
relations for creation and destruction operators with respect to a convenient reference 
state. Yet in our case the Yang-Baxter equations do not provide us with commutation 
relations, so that the algebraic Bethe Ansatz is not available. In reference pi the 



rational limit of ( |1.3|) was used to build and solve, via the algebraic Bethe Ansatz, a 
large family of isotropic multistate vertex models based on the superalgebra osp(M|2n). 

The paper is organized as follows. In Section 2, we describe the graded representa- 
tions of the TL algebra arising from orthosymplectic quantum supergroups. In Section 

3, we present the graded TL Hamiltonians with their spin correspondence. In Section 

4, the coordinate Bethe Ansatz solutions are presented with periodic, non-local and 
free boundary conditions. Finally the conclusions are reserved for section 5. 

2 Graded representations of the TL algebra 

Graded representations of the TL algebra, commuting with quantum supergroups can 
be constructed in the following way [^, ^3|. Suppose U q [XM\n] is the universal en- 
veloping superalgebra of a finite dimensional Lie superalgebra X M \ n , equipped with 
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the co-product A : U q — > U q ® U q . If now tr : U q — > End Va. is a finite dimensional 
irreducible representation with highest weight A and we assume that the decomposi- 
tion Va ® Va is multiplicity free and includes one trivial representation on Vo, then the 
projector Vq from Va ® Va onto V is a representation of the graded TL algebra. 



U k -> (g + g^ 1 )^® ... ® J® P 1<A;<JV-1 



(2.1; 



The deformation parameter g is related to Q as: 

Q + Q-^Stry^g 2 ') (2.2) 

where Str is the supertrace and p is the half-sum of the positive even roots minus that 
of the positive odd roots of X M \ n . 

Here we will consider the orthosymplectic super algebras, i.e. X M \ n = osp(M|2n). 
In order to display explicitly the Hamiltonians to be diagonalized, we will just lift the 



relevant formulas from Zhang's paper 123] 



According to the value of M the series of quantum orthosymplectic supergroups 
can be divided into four classes (see e.g. Ref. P5|): Wjosp(l|2n)], W g [osp(2|2n)], 
U q [osp(2m + l|2n)],m > 1 and W ? [osp(2m|2n)], m > 2. Namely, the vector representa- 
tions of U q [B(0,n)], U q [C(n + 1)], U q [B(m,n)} and U q [D{m,n)], respectively. 

Introducing the standard unit matrix Epq such that (Epq)rs = $pr $qs, where 
P, Q, etc. collectively denote the even and odd indices. More precisely, we have 
P = a G J with J = {1,2,...,M} or P = a e Ji with J x = {1, 2, 2n}. The 
gradation is defined as 

[P] = 0, if P € J (even) and [P] = l if P E Ji (odd). (2.3) 

and the quantum number notation is 



q x — q x 
q-q- 1 



(2.4) 



Everywhere we shall use graded-tensor product law 

(P ® Q)(\v) ® \u}) = (_l)Wl[l«>](p |„) (g, q (2.5) 
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which can be written in components as 

(P®Q)^ = (-i) m+wm P ilh Qi aj2 , 

and also the rule 

flu)® |«))t= (_i)[l«>][l«>l ( v \®( u \ 

and regard Epq as even when [P] + [Q] — (mod 2 ) and odd otherwise. 
Let us introduce a Z 2 -graded vector space with a basis 

M 

ep = (e a ,<y, a = 1,2,..., [y], a = l,2,...,ra 
where [4p] is the integer part of 4p, and a bilinear form 



(2.6) 
(2.7) 



(2.8) 



< e , e 6 >= 5 ab , a,b E J 
< ep, eQ > = < < e a , 5 a >= , a G Jo, a G J\ 
^ < 5 a ,5/3 >= -5 a p , «,/3gJi 

Then the simple roots of osp(M|2n) are particular elements of this vector space. 



(2.9) 



We extend the suffix of e P to -[f] - n < P < n + [f] by setting 



e P = 



e a 


, P = 


a <m 










, P = 


a + m, 


a < m 


a = 1, . 


.,2m 




, P = 


a < n 








. £-a = -5 a 


, P = 


a + n, 


a < n 


a = 1, . 


..,2n 



(2.10) 



Hence eo = 0. We have also introduced the index set J = Jo U J±. For each P G J, let 
— (l e a) , \8a)) £ denote the normalized weight vector having the weight ep = e a 
when P = a G Jo or ep = 8 a when P = a E J±. Thus, we can express p and a set ^4 of 
weights appearing in the vector representation n of U q [osp(M\2n)} as: 

M = 1: Vector representation of U q [B(0, n)], n>l 

Jo = {0} 



e(a) 



Ji = {±l,±2,---,±n} 
e(0) = 1 

^-l) a+n sign(a), a G Ji 

n y 

p = J2( n - a+ 2) s * 

a=l 

A={0,±6 1 ,---,±6 n } 
Q + Q' 1 = 1 - [2n\ q 



(2.11) 
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M = 2: Vector representation of U q [C(n + 1)], n > 1 

Jo = {±1} 
Ji = {±l,±2,---±ra} 
e(a) = 1, a E Jo 
e(a) = (-l) a+1 sign(a), « G Ji 

n 

p = (n — a + l)5 a — ne i 

Q=l 

.4 = {±ei,±(Ji,---,±<J n } 
Q + Q" 1 = -(?" + (2-12) 
M = 2m + 1: Vector representation of U q [B(m,n)\, m > 1, n > 1 

Jo = {0,±l,...,±m} 
Ji = {±l,±2,---±n} 
e(0) = -1, e(a) = (-l) a+m , a e J 
e( a ) = (-l) a+ri+m sign(a), a G J x 

n 1 m 1 

p = J2( n - m - a + + XX m - a + o) e " 

Q=l 2 0=1 2 

.A = {0, ±ei, ±e m , ±<Ji, • • • , ±5„} 

Q + Q- 1 = 1 - [2(n-m)], (2.13) 

M = 2m: Vector representation of U q [D(m,n)], m > 2, n > 1 

J = {±l,...,±m}, 
J 1 = {±l,±2,---±n}, 
e(a) = (-l) a , a G J 
e (a) = (-l) a+n sign(a), a G Ji 

n m— 1 

P = ^2(n - m - a + l)5 a + ^ (m - a)e a 

a=l a=l 

-4 = {±ei, ±e m _i, ±5i, • • • , ±5 n } 
Q + Q- 1 = 2 - (<f ~ m + q~ n+m )[n - m - l\ q (2.14) 

where e(P) = ±1 is a sign factor depending on the choice of X M \ n as specified 
above. 



The one-dimensional submodule of U q [osp(M\2n)} specifying the projector Vo is 
spanned by the following unnormalized singlet 

l*> = I £ e(a) q <p ' Sa> \S a ) <g> \6_ a ) + £ e(a) q <p ' €a> \e a ) <g> |e_ ) 1 (2.15) 
and its dual 

(*! = {-£ <oc) <?<'•*«> (S a \ <g> (M + £ <? <p ' ea> ® } (2.16) 



Here we have used the gradation ( |2.3| ) and the rule (|2.7[) . Using ( |2.15| ) and ( [2.16Q it is 



very simple compute the supertraces listed above 

Q + Q- 1 = q <p ' 2ta> - £ q <p ' 2Sa> - (2.17) 

Denoting the matrix unit by Epq G End V\, i.e., E ab |e c ) = 5fe c |e a ) , a,b,c £ J 
and E'q/j |<5 7 ) = <5/3 7 |5 a ), a,/3, 7 G Ji, the projector then takes the form 

1 



g + g 



-1 



£ e(a)e(/3) q <p ^ +&l3> E aP <g> E_ a _p 

a,/3eJi 

+ £ £ e(a)e(a) a^^^ ® £_ a _ a + £ aa ® £_ a _ Q ) 

0GJ1 aeJo 



+ J2 < a X b ) q <P ' ta+tb> E ab ® E„ a _ b \ . (2.18) 
a, fee Jq 



It follows from (|2.1|) that the defining relations of the TL algebra ( |1.1| ) are auto- 
matically satisfied. 

3 Graded TL quantum spin chains 

Consider a one- dimensional lattice populated with an interacting "spin " at each site 
1 < k < N. Specifically, the spin variables range over the set of weight vectors {vp \ P G 
Jo U Ji} in the following way: For M-odd, we associate a quantum spin chain of N sites 
each with spin integer s = n + m and for M-even, the spin associated is semi-integer 
s = (n+m) — 1/2. In the basis where S| is diagonal S% = diag(s, s — 1, — s + 1, — s)jt, 
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the vectors \5± a ) , a G J\ are identified with the eigenvectors of S z with eigenvalues 
±a s = ±(s — a + 1) and the vectors |e±„), |0) , a, G Jo, with the eigenvectors of 
with eigenvalues ±a s = ±(s — n — a + 1) and 0, respectively. Note that this map recasts 
even and odd suffixes into the eigenvalues of S z : 

a s G j{ s \ {a s , 0} G J (s) and J (s) = J (s) U J x (s) = {s, s - 1, -s + 1, -s}. (3.1) 

Thus the Hilbert space is an iV-fold tensor product V\ ® N . The Hamiltonians 
associated with the representations of the graded TL algebra are sums of the TL 
operators: 

N-l 

+ b.t. (3.2) 

k=\ 

where Uk = Uk,k+i operates in a direct product of Va at positions k and k + 1. In 
general the boundary terms b.t. break translational invariance, reflecting the non- 
cocommutativity of the co-product. So, very special boundary terms must be consid- 
ered when we seek quantum (super)group invariance of H. In particular, one possi- 
bility is to consider free boundary conditions, i.e., b.t. = 0. In the next Section we 
consider another boundary term, first presented in the framework of the coordinate 
Bethe Ansatz in |[26|| , here named Martin's boundary conditions \2%, for which H is 



quantum supergroup invariant. 

Besides the explicit matrix elements, available from ( |2.11| )-( f2.14| ) and ( |2.18 ), the 



local Hamiltonians densities are expressible in a more convenient form using ( |2.15| ) and 

u k = mm) k (3.3) 

For instance, let us consider W g [osp(l|2n)] Hamiltonians. In particular, for the simplest 
case, i.e., the W g [osp(l|2)] case, the spin correspondence gives us a spin-1 Hamiltonian. 
From (|2~TT| ), (|2~TB| ) and (fO§, we have 



m = g- 1/2 |^i)®|5-i)-g 1/2 |5-i)®|5i) + |0)®|0) 

= -g-^^l^^l+gVa^l^^l + ^l®^! (3.4) 
In the basis where Sf. is diagonal with eigenvectors |+, k) , |0, k) , | — , k) and eigenvalues 
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+1, 0, —1, Uk is the 9 by 9 matrix acting on \a, k) <S> \b, k + 1), a, b = +, 0, 

/ \ 

u k = mm k , k+1 = I o v o I (3.5) 



k,k+l 



where is the 3 by 3 zero matrix and V is 



q 1 q 



(+,-1 

V= (0,0| I -q- 1 ' 2 1 
(-,+1 I 1 



1/2 j 



gl/2 



(3.6) 



In general, for Wjosp(l|2n)] the correspondent spin is s = n and the local Hamil- 
tonian U k is the (2n + l) 2 by (2n + l) 2 matrix acting on \a,k) (g) |6, k + 1), a,b = 



n, 1,0, -1, -n: 



17* 



/ 











V 





o \ 







(3.7) 



Vo ... o o o ■•• o; kk+1 

where is the (2n + 1) by (2n + 1) zero matrix and V is 



V 



-q)- n Q n A nl Q n \ 

e n {- q ) n A nl (- q ) n e n J 



(3i 



in the basis {|n, — n) , |1, — 1) , |0, 0) , |— 1, 1) , |— n, n)}. Here 6 n is the n by n 
matrix 

(©n) y = {-q)- n - l+l+ \ i,j = 1,2, -,n (3.9) 
and A n i is the n by 1 matrix 



(Au)a 



{-l) n - 1 g - n+i " 1 /2 j i = l,2,...,n. 



(3.10) 



A^j stands for the transposed of A nl . 

Finally in view of the grading the basis vectors of the module V\ can be identified 
with the electronic states (see e.g. Ref. [IB|). Thus, one can get alternative expressions 
for the local Hamiltonians U k in terms of the canonical fermion operators. These 
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expressions are not presented here since they are awkward and do not help us in the 
diagonalization of H. 

These Hamiltonians appear to be new although due to the TL equivalence , they 
are expected to possess the same thermodynamics as the ZY g [osp(l|2)] spin-1 chain with 
appropriate coupling. Thus, we also expect that the W 3 [osp(l|2)] spin-1 model plays 



the same role that the spin-1/2 XXZ chain plays in non-graded case [21 



Having now built common ground for all W g [osp(M|2n)] Hamiltonians, we may 



follow the steps of [18| to find their spectra 



4 The coordinate Bethe Ansatz 

Using a spin language one can say that the graded TL Hamiltonians are spin-zero 
projectors. They actually are the projectors on total two-site spin zero. 

Since H commutes with the total spin operator , its spectrum can be classified 
in sectors which are defined by the eigenvalues of the operator number r = sN — Sy, 
where the total spin operator is defined by 

N 

S^ = J2l® (k ~ 1] ®S z k ®I® {N - k) . (4.1) 

k=l 

Therefore, there exists a reference state ^q, satisfying H^/ = E ^ , with E = 0. 
We take \&o to be = \s s s ■ ■ ■ s). It is the only eigenstate in the sector r = and 
all other energies will be measured relative to this state. 

In every sector r there are eigenstates degenerate with \Ev They contain a set of 
impurities. We call impurity state obtained by lowering some of the \s, k)'s, such that 
the sum of any two neighboring spins is non-zero. Since if is a projector on spin zero, 
all these states are annihilated by H. In particular, they do not move under the action 
of H, which is the reason for their name |28|| . 

Nothing interesting happens in sectors with r < 2s. In sector r = 2s, we encounter 
the situation where the states \—x,k) and \x,k±l), x G J^ s \ occur in neighbor- 
ing pairs. They move under the action of H, i.e., the sector r = 2s contains one 
free pseudoparticle. In general, for a sector r we may have p pseudoparticles and 
N s -.i, iV s _ 2 , iV_ s+ i impurities of the type s — l,s — 2, s + 1, respectively, such 
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that 

2s-l 



2sp + J2 xN *-x- (4-2) 



x=l 



The main result of this section is to show that H can be diagonalized in a convenient 
basis, constructed from products of single pseudoparticle wavef unctions. The energy 
eigenvalues will be parametrized as a sum of single pseudoparticle contributions. 

In sector r = 2s, the states \k{—(3,(3) >= \s s- • -s -[3(3 s - ■ - s > span the corre- 

k 

sponding eigenspace of H, where k = 1,2, JV — 1 and (3 E J^ s \ We seek eigenstates 
of H which are linear combinations of these vectors. It is very convenient to consider 
the vector 



\m) = E <l)e(a s )q< 5 ^>\k(-a s ,a s )) 

+ e(l)e(a s )q <S ^>\k(-a s ,a s )). (4.3) 



Here we observe that this linear combination corresponds exactly to the first row vector 
of matrix Uy.. For instance, we can recall (|3.6|) to write down | for osp(l|2): 

\VL{k)) = ) + +:-+ + •••+> + q- 1 ' 2 \ + + : 00+ ••• + > 

' fc k (A A) 

- q" 1 ] + + +- + •••+ >, {k = 1,2,..., N -1) 

fc+i 

The choice of \Q(k)) is very special. It is an eigenvector of Uf. with eigenvalue equal 
to supertrace (p.2[) 

U k \tt(k)) = (Q + Q- 1 )\n(k)), (4.5) 
and the action of Uk±i on is easily done: 

E7 fc _! |n(A:)) = \Q(k - 1)) \Q(k)) = \n(k + 1)) 

U k \Q(m)) = k^{m±l,m} 



(4.6) 



It should be emphasized that although the linear combination ( [4.3| ) is different 
for each model, the action of Uf. is always given by fl4.5| ) and ( |4.6| ). Therefore, all 
orthosymplectic Hamiltonians (|3.2|) can be diagonalized in a similar way affording a 
considerable simplification in their diagonalizations when we compare with the calculus 



used in the usual spin basis [28 . 
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4.1 Periodic boundary conditions 
4.1.1 One-pseudoparticle eigenstates 

We will now start to diagonalize H in every sector. Let us consider one free pseudopar- 
ticle as a highest weight state which lies in the sector r = 2s 

* 2s = J2A(k)\n(k)). (4.7) 



where \Q(k)} is given by ( (Of ). Using the eigenvalue equation H ty 2s = E 2s ty 2s} one 
can derive a complete set of equations for the wavefunctions A(k). 

When the bulk of H acts on \Q(k)) it sees the reference configuration, except in the 
vicinity of k where we use (|4.5| ) and (|4.6|) to get 

H\a(k)) = (Q + Q~ l )\n(k)) + \n(k-i)) + \n(k + i)) 

2 < £;<iV-2 (4.8) 



Substituting ( |4.8|) in the eigenvalue equation, we have 

(E 2s - Q - Q- l )A(k) = A(k - 1) + A{k + 1) 

2 < k<N-2 (4.9) 

Here we will treat periodic boundary conditions . They demand t.b = Z7jv,i, implying 
A(k + N) = A(k). This permits us to complete the set of equations (|4.9|) for A(k) by 
including the equations for k = 1 and k — N — 1. Now we parametrize A(k) by plane 
wave A(k) = A£ k to get the energy of one free pseudoparticle as: 

£ N = 1 (4.10) 

Here £ = e* e , 9 being the momenta determined from the periodic boundary to be 
9 = 2ttI/N, with / an integer. 

4.1.2 One-pseudoparticle and impurities 

Let us consider the state with one pseudoparticle and one impurity of type (s — 1), 
which lies in the sector r = 2s + 1. We seek eigenstates in the form 
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*Wi(6,6) = E fe) |ni(A:i, fc 2 )> + A 2 (A;i, A; 2 ) |0 2 (A:i, A; a )>} (4.11) 

fci<fc 2 

We try to build these eigenstates out of translationally invariant products of one pseu- 
doparticle excitation with parameter £2 and one impurity with parameter £1: 

*Wi(£i,6) = IO-i)(£i)> x <M6) + * 2s (6) x |(s-i)(eO> (4-12) 



Using one-pseudopaticle eigenstate solution (|4.3|) and comparing this with ( [4.11 

get 



we 



\tti{h,k 2 )) = E e{l)e{a s )q <5l+5as ' p> \k 1 {s-l),k 2 {-a s ,a s )) 

+ E q <5l+ea °' p> Ms-l),k 2 (-a s ,a s )) (4.13) 



and similar expression for \Q 2 (ki, k 2 )), for which the impurity (s — 1) lies at position 
k 2 . Moreover, the wavefunctions are parametrized by plane wave 

A 1 (k 1 , k 2 ) = A x && , A 2 {kx, k 2 ) = A 2 &£\ (4-14) 

Periodic boundary conditions A x (k 2 , N + hi) = A 2 (ki, k 2 ) and A^N + ki, N + k 2 ) = 
Ai(ki, k 2 ), i — 1,2 imply that 

A^ = A 2 , ^ = (a/ = l (4.15) 

When H now acts on ^ 2s+ i, we will get a set of coupled equations for A^ki, k 2 ), 
i = 1,2. We split the equations into far equations, when the pseudoparticle does not 
meet the impurity and near equations, containing terms when they are neighbors. 

Since the impurity is annihilated by H, the action of H on ( |4. 1 1[ ) in the case far 
(i.e., k 2 > ki + 2), can be written down directly from 



(E 2s+1 -Q- Q- 1 ) A x {k x , k 2 ) = A 1 (k 1 , k 2 - 1) + A 1 (k 1 , k 2 + 1) (4.16) 

and a similar set of equations with A 2 (the pseudoparticle at position k x ). Using the 
parametrization ( |4.14|) , these equations will give us the energy eigenvalues 



E 2s+1 = Q + Q- 1 + ^ 2 + ^ 2 1 (4.17) 
12 



To find £2 we must consider the near equations. First, we compute the action of H on 
the coupled near states \Qi(k,k + 1)) and \Q 2 (k, k + 2)): 



H |fii(A;, k + 1)) = (Q + Q-^lfiifofc + l)) 

+ |n 1 (fc,fc + 2)) + |fi 2 (fc,fc + 2)> (4.18) 

and a similar set of equation for \fl 2 (k, k + 2)). The last terms in these equations tell 
us that a pseudoparticle can propagate past the isolated impurity, but in so doing 
causes a shift in its position by two lattice site. Substituting (|4.18j) into the eigenvalue 
equation, we get 

(E 2s+1 -Q- Q- 1 ) A^k, fc + 1) = Ai(k, k + 2) + A 2 {k, k + 2) (4.19) 
These equations, which are not automatically satisfied by the ansatz ( [4. 14 ), are equiv- 



alent to the conditions 

A t (k,k) = A 2 (k,k + 2) . (4.20) 



obtained by subtracting Eq.( |4.19i ) from Eq. ( |4. 16|) for k\ = k , k 2 = k + 1. The condition 



( |4.20j ) requires a modification of the amplitude relation (|4.15|) : 

f x = g = ^ 2 => = 1 or = 1 (4.21) 

Putting & = e l9i , i — 1,2, it means cos(A^ — 2)9 2 = cos26>i. Hence 

2-nm±A'nl/N , tA , 

2 = — — , I and m integers. (4.22) 

In other words, \?2s+i(6b £2) are eigenstates of H with energy eigenvalues given by 
E 2s+ \ = Q + Q~ x + 2 cos 9 2 . Note that when Nm it 21 is a multiple of (N — 2) we get 
states which are degenerate with the one-pseudoparticle states ^ 2s , which lie in the 
sector r = 2s. 

In the sectors r = 2s + I, I > 1 we also will find states, which consist of one 
pseudoparticle with parameter interacting with I impurities, distributing according 



to (4^), with parameters £j, i — 1, 2..., /. 
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The energy of these states is parametrized as in Q4.17| ) and 6+1 satisfies the condition 
( Egg ) with f = & • • • £j 6+1- Jt involves only 6+1 and £ imp = 6 6 • • • 6, bein g therefore 
highly degenerate, i.e. 

fei £2 •••£? = ! (4-23) 

This is to be expected due to the irrelevance of the relative distances, up to jumps 
of two positions via exchange with a pseudoparticle. Moreover, these results do not 
depend on impurity type. 

4.1.3 Two-pseudoparticle eigenstates 

The sector r = 4s contains, in addition to the cases discussed above, states which 
consist of two interacting pseudoparticles. We seek eigenstates in the form 

¥4,(6,6)= E ^1,^)10(^^2)) (4.24) 
fci+i<fc 2 

Again, we try to build two-pseudoparticle eigenstates out of translationally invariant 
products of one-pseudoparticle excitations at k± and k 2 {k 2 > k\ +2) : 

¥4,(6,6) = *2,(6) x ¥ 2 ,(6) + ¥ 2s (6) x ¥ 2s (6) (4.25) 
Using one-pseudoparticle solution ( f4.7| ) and comparing ( |4.24j ) with (|4.25|) , we get 



A(k u k 2 ) = A u && + A^il 2 (4-26) 



and 



|ni(A:i,*:2)> = 

= E e(a s )e(f3 s )q <2Sl+Saa+&Ps ' P> \h(-ot s , a s ), k 2 (-(3 s , &)> 



E e(aMPs)q <25l+5 ^' p> \h(-a s , a s ), k 2 (-b s , b s )) 
E e(a s )e(P s )q <25 ^ +5 ^ > \k 1 (-a s ,a s ),k 2 (-P s ,P s )) 
+ E <aMbs)q <2Sl+ea ° +eb °' p> \h(-a s ,a s ),k 2 (-b s ,b s )) (4.27) 



a s e4 s) ,/3 a eJ[ s) 



for k 2 > fci +3. 
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Periodic boundary conditions A{k 2l N + ki) = A(ki, k 2 ) and A(N + ki,N + k 2 ) = 
A(k 1 ,k 2 ) imply 

Ai2$ = A 21 and £"=1 (4.28) 

where £ = £1^2 (£» = e * e % * = 1,2) and the total momentum is Q\ + #2 = 2irl/N, with Z 
integer. 

Applying i7 to the state of (|4.24j) , we obtain a set of equations for the wavefunctions 
A(ki,k 2 ). When the two pseudoparticles are separated, (k 2 > k\ + 3) these are the 
following /or equations: 

(E 4s -2g-2g- 1 )A(A; 1 ,A;2) = A(h - l,k 2 ) + A(h + l,k 2 ) 

(4.29) 

+ ^(/^/ca-lJ+^/^/ca + l) 

We already know them to be satisfied, if we parametrize A{k\,k 2 ) by plane waves 
( |4.26| ). The corresponding energy eigenvalue is 

E 4s = 2Q + 2Q- 1 + 6 + + 6 + (4-30) 

The real problem arises of course, when pseudoparticles are neighbors, so that they 
interact and we have no guarantee that the total energy is sum of single pseudoparticle 
energies. 

Acting of H on the near states gives the following set of equations 
H\Q(k,k + 2)) = 2(Q + Q- 1 )\n(k,k + 2)) + \Q(k-l,k + 2)) 

(4.31) 

+ \Q(k, k + 3)) + U k+1 \Q(k,k + 2)} 

Before we substitute this result into the eigenvalue equation, we observe that some 
new states are appearing. In order to incorporate these new states in the eigenvalue 
problem, we define 

U k+1 \Q(k, k + 2)) = \Q(k, k + l)) + \Q(k + l,k + 2)) (4.32) 

Here we underline that we are using the same notation for these new states. Applying 
H to them we obtain 

H\Q(k,k + l)) = (Q + Q- 1 )|fi(fc,fc + l)) + |fi(/c-l,/c + l)) 

(4.33) 

+ \Q(k,k + 2)} 
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Now, we extend ( [4.24| ), the definition of , to 



*4.(6,&)= E A^k^Mh,^) (4.34) 
fci<fc 2 



Substituting (|4.31|) and ( |4.33| ) into the eigenvalue equation, we obtain the following set 
of near equations 

(E is -Q- Q- 1 ) A(k, k + 1) = A(k - 1, k + 1) + A(k, k + 2) (4.35) 



Using the same parametrization ( |4.26| ) for these new wavef unctions, the equation (|4.35|) 



gives us the phase shift produced by the interchange of the two interacting pseudopar- 
tides 

Ai 2 " l+^ + (Q + Q-l)^ 1 • > 

We thus arrive to the Bethe Ansatz equations which fix the values of £i and £2 in the 
energy equation ( |4.30| ) 

^ 2 i + £ + (Q + Q-% 

= (U2) N = l (4.37) 
4.1.4 Two-pseudoparticles and impurities 

In the sectors r > 4s, in addition the cases already discussed, we find states with two 
interacting particles and impurities. Let us now consider states with two pseudoparti- 
cles and one impurity, for instance, of type s — 1. Theses eigenstates lie in the sector 
r = 4s + 1 and we seek them in the form 

#4 S +i(£i,£ 2 ,£3) = E Mkuhik^in^k^h)) 

k 1 +l<k 2 <k 3 ~2 

+ e A 2 (h, k 2 , h) \Q 2 (h, k 2 , h)) 
fci+i<fc 2 <fc 3 

+ E A 3 {h,k2,h)\n 3 {h,k2,h)) (4.38) 

fc 1 +l<fe 2 <A:3-l 

In Aj(ki, &2, ^3) the index j = 1, 2, 3 characterizes the impurity position. 

Comparing (|4.38|) with the states build from the translationally invariant products 
of two-pseudoparticles with parameters £2 and £3 and one-impurity with parameter £1: 

*4, + i(£i,6,&) = |(s-i)(£i))x* 2s (6)x* 2s (£ 3 ) 

16 



+ |(s-l)(£i))x* 2s (£ 3 )x* 2s (£ 2 ) 

+*2.(es)x|( a -i)(ei)>x* 2 .(6) 

+ ^2 S (6)x^ 2s (e 3 )x|( S -l)(a)) 

+#2,(6) x* 2a (6) x |(s - l)(Ci)> (4-39) 

we get 

fc 2 , fa)) = 

= Yl e(a s )e(p s )q <25 ^ +s ^ > \fa(s-l), k 2 (-a s , a s ), fa(-(3 s , &)> 

£ e(a s )e(Ps)q <25l+5a ° +eb °' p> \fa(s-l), k 2 (-a s , a,), fa(-b„ b s )) 
]T e(a s )e(^)g <2<5l+e - + ^^ - 1), fc 2 (-a a , a s ), fc 3 (-ft, &)> 

a s e4 s) ,/3 sG J< s) 

+ Yl e(a s )e(bs)q <25l+ea ° +eb °> p> \fa(s - 1), k 2 (-a s , a,), fa(-b„ b,)) 

(s) 

(4.40) 

and similar expressions for \Q 2 (fa, k 2 , fa)) and \Q 3 (ki, k 2 , fa)). This comparing also 
gives us the parametrization for the wavef unctions A^(k±, k 2 , fa) : 

Ai{faMM = A 23 ststst + 

A 2 (fa, k 2 , fa) = A 213 &&& + A 231 &&& 

A s (fa,k 2 ,fa) = A 312 e l s e 2 i &+A S21 e l 3 e 2 2 e 3 i - (4.41) 

Periodic boundary conditions read now 

A } (N + fa, N + k 2 , N + fa) = A s {fa,k 2 ,fa), 



A s (k 2 ,fa,N + fa) = A j+1 (fa,k 2 ,fa) 

j = 1,2,3 (mod3) (4.42) 

which imply that 

p N A 123 _ A 132 N _ A 3 i2 _ A 2 13 

SI — ~1 — ~1 J S2 _ ~~A — ~~A ) 

^312 ^321 ^231 ^132 
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rN 
S3 



.4 



321 



Ai 



^213 ^123 

Action of if on the state ^4 S +i gives the following set of far equations: 



(4.43) 



(E 4s+1 -2Q-2Q- 1 ) A 1 {k 1 ,k 2 ,k 3 ) = A 1 {k 1 ,k 2 -l,k 3 ) + A 1 {k 1 ,k 2 + l,k 3 ) 

+A 1 (k 1 , k 2 , k 3 -l) + A 1 {k 1 , k 2 , k 3 + 1) 

(4.44) 

and a similar set of eigenvalue equations for A 2 (ki,k 2 ,k 3 ) and A 3 (kx, k 2 , k 3 ). The 



parametrization ( |4.41| ) solves these far equations provided: 

E, s+l = 2Q + 2Q- 1 + 6 + + 6 + 4V 1 



(4.45) 



Taking into account the near equations we must split them in three different neigh- 
borhood: (i) impurity neighbors of separated pseudoparticles, (ii) impurity far from 
neighbors pseudoparticles and (iii) when impurity and pseudoparticles share the same 
neighborhood. 

In the case (z) we consider the second pseudoparticle far and follow the steps for 
the case of one-pseudoparticle with impurity eigenstates presented above. Thus, the 
near equations can be read off from ( f4.35| ) 



(E 4s+1 -2Q -2Q- 1 )A 1 (k,k + l,k 3 ) = A 1 (k,k + l,k 3 -l)+A 1 (k,k + l,k 3 + l) 

+A 1 (k, k + 2, k 3 ) + A 2 (k, k + 2, k 3 ) 

(4.46) 

and a similar set of equations coupling A 2 and A%. It follows from the consistency 
between ((03) and (KW) that 



A 1 (k, k, k 3 ) = A 2 (k,k + 2,k 3 



(4.47) 



and similar identification between A 2 and A s . The plane waves fl4.41|) solve these 
identifications provided: 

^4l23 ^132 ^231 ^213 



A 2 i 3 A 23 i A 32 i A su 
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(4.48) 



For the case (ii) we can derive the near equations from those of two-pseudoparticles 
case. Keeping the impurity far and following the steps ( |4.29| )-( p.36| ), we get 



(E As+1 -Q- Q- 1 ) k, k + 1) = A 1 (k 1 , k - 1, k + 1) + A x {k x , k, k + 2) (4.49) 



and a similar set of equations for A 2 and v4 3 . The case (iii) is obtained from ( 4.49|) for 
kx = k-l. 

The parametrization ( 4.41]) solves these cases provided: 



A123 = _ i + &6 + (Q + <r% f4 5 ) 

Matching the constraint equations ( |4.50| ), ([4.48 ) and ( |4.43| ) we arrive to the Bethe 



equations 

gtl = - \Vf^(ntn~~Z a > «^ & = 2 >3 (4.51) 
^ = (666)" = i, Cf ~ 4 = 1- (4.52) 

The origin of the exponent (AT — 4) in the impurity parameter can be understood by 
saying that after the two pseudoparticles propagate past the impurity, the position of 
impurity is shifted by four lattice sites. 

Next, we can also find eigenstates with two pseudoparticles and more than one 
impurities. They can be described in the following way: Let us consider an eigen- 
state with I > 1 impurities with parameters £i, £ 2 , ■ • ■ > 6 an d two pseudoparticles with 
parameters and £z + 2- The energy eigenvalue is 

E r = 2Q + 2Q- 1 + + C+i + 6+2 + 6+2 (4-53) 

and the Bethe equations 



f JV / r 2/ r 2 f 2 
SZ+1S1S2 ' ' ' Q 



i + 6+16+2 + (g + g- 1 )^+i 



1 + 6+16+2 + (Q + Q- 1 )6+ 2 

ef- 4 = 1, a = 1,2,...,/ (4.54) 



Moreover, £ w = 1 with f = • • -66+16 



Z+2 
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(4.56) 



4.1.5 Three-pseudoparticle eigenstates 

In the sector r = 6s, in addition to the previously discussed eigenstates of one and 
two pseudoparticles with impurities, one can find eigenstates with three interacting 
pseudoparticles with parameters £1,^2 and £3. We start seek them in the form 

* 6s (£i,6,6)= E v4(/ci, k 2 , k 3 ) k 2 , k 3 )) (4.55) 

fci+2<fc 2 <fc3-2 

where \fl(ki, k 2 , k 3 )) = <S>f =1 \Q(ki)) . The corresponding wavefunctions 
A(k u k 2 , k 3 ) = A 123 ei^2 2 & 3 + A^ 1 ^^ 2 + A 213 ei 2 && 

+^23l£l <52 <53 + A 312<Sl <?2 <53 + A 321^1 <?2 ?3 

satisfy the periodic boundary conditions 

A(k 2 , k 3 , N + ki) = A(k u k 2 , k 3 ), A(N + ki,N + k 2 , N + k 3 ) = A(k u k 2 , k 3 ) (4.57) 
which imply that 

^123 _ ^132 _ ^312 _ ^213 

1 ^312 ^321 ' ~ ^231 ^132 ' 

g = ^1 = ^1, ^ = (^ 3 ) N = l (4.58) 

^213 ^123 

These relations show us that the interchange of two-pseudoparticles is independent of 
the position of the third particle. 

Applying H to ( f4.55|) , we obtain a set of equations for the wavefunctions A(ki, k 2 , k 3 ). 
When the three pseudoparticles are separated, (ki + 2 < k 2 < k 3 — 2), we get the fol- 
lowing far equations: 

(Ees-ZQ- SQ-^Aih^^h) = A{h - 1, k 2 , k 3 ) + A{h + 1, k 2 , k 3 ) 

+A(k u k 2 - 1, k 3 ) + A(ki, k 2 + 1, k 3 ) 
+A(k u k 2 , k 3 - 1) + A(ki, k 2 , k 3 + 1) 

(4.59) 

It is simple verify that the wavefunctions ( [4.56| ) satisfy these far equations provided 



Ess = Y.{Q + Q~ 1 + tn + C 1 } (4-60) 



n=l 
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Applying H on the near states we get the following set equations: 

H\n(k 1 ,k 1 + 2,k 3 )) = (2Q + 2Q- 1 )\Q(k 1 , k 1 + 2,k 3 )) + 1^(^-1, h + 2,k 3 )) 

+ \Q(k u h + 3, fc 3 )> + \Q(k u k x + 2,k 3 - 1)> 
+ \Q{k u ki + 2,k 3 + 1)) + U kl+ i \Q{h, h + 2, k 3 )) 

(4.61) 

for k 3 > ki + 4, which correspond to the meeting of two pseudoparticles at the left of 
the third pseudoparticle, which is far from of the meeting position. 

H\n{k u ki,k 2 + 2)) = (2Q + 2Q- 1 )\n(k 1 ,k 2 ,k 2 + 2)} + \n(k 1 -l,k 2 ,k 2 + 2)} 

+ + 1, k 2 , k 2 + 2)) + \n(k h k 2 -l,k 2 + 2)) 
+ \Q{k u k 2 , k 2 + 3)) + U k2+1 \Q(h, k 2 , k 2 + 2)) 

(4.62) 

for k 2 > ki + 2, which correspond to the meeting of two pseudoparticles at the right 
of the far pseudoparticle. Moreover, there is one set of equations which correspond to 
the meeting of three pseudoparticles 

H\tt(k,k + 2,k + 4)} = (Q + Q~ l ) k + 2, k + 4)) + \£t(k - l,k + 2, k + 4)) 

+ \Q(k, k + 2,k + 5)} + U k+ i \Q(k, k + 2, k + 4)) 
+U k+3 \n(k,k + 2,k + A)} (4.63) 

In deriving these equations new states made their debut. In order to incorporate these 
new states in the eigenvalue problem we define: 

E/ fcl+ i|fl(A;i,A; 1 + 2,A;g)) = \Q{k u h + 1, k 3 )) + \Q{h + 1, h + 2, k 3 )) 

u k2+x \n{k^k 2 ,k 2 + ^)) = \n(k 1 ,k 2 + i,k 2 + 2)) + \n(k 1 ,k 2 ,k 2 + i)) 

U k+1 \Q(k,k + 2,k + 4:)) = \Q(k, k + 1, k + 4)) + \Q(k + 1, k + 2, k + 4)) 
U k+3 \tt(k,k + 2,k + 4:)) = \n(k,k + 3,k + A)) + \n(k,k + 2,k + 3)) (4.64) 

Applying H to these new states the result can be incorporated to the eigenvalue 
problem provided the definition of ^q s Q4.55 ) is extended to 



*6,(£i,6,&)= E A{k u k 2 MP{Kk 2 ,k 3 )) (4-65) 
fci<fc 2 <fc 3 
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After this we are left with three meeting equations 

(E 6s -2Q- 2Q- 1 )A(k 1 , h + 1, k 3 ) = A(kx - 1, A* + 1, k 3 ) + A(h, k x + 2, k 3 ) 

+ A(k h h + 1, k 3 - 1) + A(kx, h + 1, k 3 + 1) 

(4.66) 

for &; 3 > &i + 2, 

(£ 6s - 2Q - 2Q- 1 )A(k 1 , k 2 , k 2 + 1) = A; 2 — 1, fc 2 + 1) + A(fci, fc 2 , A; 2 + 2) 

+ A(k t - 1, fc 2) k 2 + 1) + + 1, fc 2 , fc 2 + 1) 

(4.67) 

for ki + 2 < k 2 and 

(E 6s -Q- Q~ 1 )A{k, k + 1, k + 2) = A(k - 1, A; + 1, A; + 2) + A(Jfe, fc + 1, A; + 3) 



It is easy to verify that the parametrization ( |4.56| ) and ( |4.60| ) solve these equations 
provided 



^123 


^231 


1+62 


+ (Q + Q~ 




^213 


^321 


I+62 


+ (Q + Q~ 




^4l32 


^213 


1+^13 


+ {Q + Q~ 




^231 


^312 


1+^13 


+ {Q + Q~ 




^312 


^4l23 


1+63 


+ {Q + Q~ 




^321 


^4l32 


1+63 


+ (Q + Q~ 





(4.68) 

where 6j — £i£r Matching these constraints and the periodic boundary conditions 
Q4.58Q we get the Bethe Ansatz equations 



3 



s = n -y 



a = 1,2,3 



fates)" = 1 (4.69) 
4.1.6 General eigenstates 

The generalization to any r is now immediate. Since the Yang-Baxter equations are 
satisfied, there is only two-pseudoparticle scattering (using the ^-matrix language). 
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Therefore, neighbor equations, where more then two pseudoparticles become neighbors, 
are nor expected to give any new restrictions. For instance, in the sector r = 6s, we saw 
that the interchange of two-pseudoparticles is independent of the position of the third 
particle. Thus, in a sector with p pseudoparticles we expect that the p-pseudoparticle 
phase shift will be a sum of p ^ p ~ 1 - ) two-pseudoparticle phase shift. The energy is given by 
the sum of single pseudoparticle energies. The corresponding Bethe Ansatz equations 
depend on the phase shift of two pseudoparticles and on the number of impurity. For 
a generic sector one can verify that no different neighborhood those discussed above 
can appear. So, no additional meeting conditions will be encountered. Thus, we can 
extend the previous results to the p -pseudoparticle states in the following way: In a 
generic sector r with / impurities parametrized by 662 • • • 6 and p pseudoparticles with 
parameters 6+16+2 ■ • • 6+ P , the energy is 

E r = £ {Q + Q-'+tn + tn 1 } (4.70) 
n=l+l 

with £ n determined by the Bethe ansatz equations 

,iV,2,2 ,2 rf f 1+ tea + (<? + <?-% ! 

^ J-i+il i+ a* +«?+<?-% / 

e~ 2p = i, c =i,2,...,/ 

Z N = 1, £ = 66 •• -66+i6+2 •••£+*>■ (4.71) 

The energy eigenvalues and the Bethe equations depend on the deformation parameter 
q, through the relation (|2.2|): 



r 1 - [2n) q for B(0,n) 

O + o- 1 - -(?" + 9- n )[n-l] 9 for C(n+1) , , 

^ + ^ | l-[2(n-m)] g for B(m,n) K > 

[ 2 - (q n - m + q' n + m ) [ n - m - l) q for D(m, n) 

We obtained thus the spectra with periodic boundary conditions of quantum spin- 
chain models, arising as representations of the graded Temperley-Lieb algebra. As 
expected, all these models have equivalent spectra up to degeneracies of their eigen- 
values. For a suitable sorting of the parameters 6, one can insure that the spectra of 
lower-r sectors are contained entirely in the higher-r sectors. 



23 



4.2 Martin's boundary conditions 

It is the purpose of this subsection to present and solve, via coordinate Bethe ansatz, 
the quantum supergroup invariant closed TL Hamiltonians which can be written as 



N-l 

H=Y.U k +U Q (4.73) 
fe=i 

where Uk is here a graded TL operator presented in the Section 2, and t.b. = IAq is 
non-local term defined through of a operator G which plays the role of the translation 
operator 

U^GUn^G- 1 , G=(Q-U 1 )(Q-U 2 )---(Q-U N - 1 ) (4.74) 

satisfying [H, G] = and additionally invariance with respect to the quantum super- 
algebra. The operator G shifts the Uk by one unit GU^G -1 = U^+i and maps U into 
Ui, which manifest the translational invariance of H. In this sense the Hamiltonian 
( |4.73| ) is periodic. 



The action of the operator G on the states \Q{k)) can be easily computed using 



75J) and (|4.6|): It is simple on the bulk and at the left boundary 



G \Q(k)) = -Q N - 2 \n(k + 1)) , l<£;<iV-2 (4.75) 
but manifests its nonlocality at the right boundary 

N-l 

G \Q(N - 1)) = Q N - 2 £ {-QY k \n{N - k)) (4.76) 

k=l 

Similarly, the action of the operator G^ 1 = (Q^ 1 — Un-i) • • • (Q^ 1 — U\) is simple on 
the bulk and at the right boundary 

G- 1 \Q(k)) = -Q- N+2 \Q(k - 1)) , 2<A;<iV-l (4.77) 
and non-local at the left boundary 

N-l 

G- 1 |fi(l)> = Q- N+2 ^ (~Q) k \m) ■ (4-78) 

k=l 

Now we proceed the diagonalization of H as we did for the periodic case. 
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4.2.1 One-pseudoparticle eigenstates 

Let us consider one free pseudoparticle which lies in the sector r = 2s 

N-l 

* 2s = E A (k) TO) . (4.79) 
The action of the operator U. = Y.k=i Uk on the states \£l{k)) can be computed using 



(H): 



zy|o(i)) = (g + g- 1 )|fi(i)) + |^(2)) 

U\n(k)) = (g + g- 1 )|fi(A;)) + |fi(A;-l)) + |fi(A; + l)) 

for 2 < k < N — 2 

U\Q{N-1)) = (Q + Q- 1 )\Q(N-l)) + \Q(N-2)). (4.80) 



and using (|4.75|) - (pL78|) one can see that the action of Uq = GU^^iG 1 vanishes on the 
bulk 

U \Q(k))=0 , 2<k<N-2 (4.81) 
and is nonlocal at the boundaries 

N-l iV-1 

Wo |^(i)) = - E (~Q) k \m) , Wo [o(at - 1)) = - E (-g)- w+fc |n(fc)) . (4.82) 

k=l k=l 

which are connected by 

U \Q(N - 1)) = (-Q)~ N U |fi(l)> . (4.83) 

From these equations we can understand the role of Uq: Although the Hamiltonian 
( [4.73j ) is a global operator, it manifests the property of essential locality. From the 
physical point of view, this type of models exhibit behavior similar to closed chains 
with twisted boundary conditions. 

Before we substitute these results into the eigenvalue equation, we will define two 
new states 

\n(0)) = U |0(1)> , |0(iV)) = U \Q(N - l)) (4.84) 
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to include the cases k = and k = N into the definition of ty 2s , equation ( |4.79| ). 
Finally, the action of H = IA + Uq on the states \fl(k)) is 



H\n(o)) = (Q + q- 1 ) \n(o)) + (-Qf \n(N - 1)) + 
H\n(k)) = (Q + q- 1 ) \n(k)) + \n(k - 1)) + \n(k + 1)) 

for 1 < k < N — 2 

h\q(n-i)) = (Q + q- 1 ) \n(N - 1)} + \n(N - 2)} + (-Q)- N \n(o)} 

H\Q(N)) = (Q + Q- 1 ) \n(N)} + \n(N - 1)} + (-Q)- N (4.85) 

Substituting these results into the eigenvalue equation H^/ 2s = E 2s \&2s we get a 
complete set of eigenvalue equations for the wavefunctions 



E 2s A(k) = (Q + Q- v )A{k) + A(k - 1) + A(k + 1) 

fori < k < N - 1 (4.86) 

provided the following boundary conditions 

(-Q) N A(k) = A(N + k) (4.87) 

are satisfied. 

The plane wave parametrization A(k) = A^ k solves these eigenvalue equations and 
the boundary conditions provided: 

E 2s = Q + Q-' + Z + C 1 

t N = (~Q) N (4.88) 

where ^ = e ie and 8 being the momentum. Therefore the coordinate Bethe Ansatz 
technique is the same used in the previous section for the periodic case. Differences 
arise from the boundary terms. Thus, let us keep the bulk's results of the previous 
section and only present the results due to the action of U 
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4.2.2 Two-pseudoparticle eigenstates 

The action of Uq does not depend on the pseudoparticles are neither separated nor 
neighbors. It vanishes in the bulk 

|n(ifei, A: 2 )> = for h ^ 1 and k 2 ^ N- 1, (4.89) 

and is different of zero at the boundaries: 

Z^|n(l,Ar 2 )> = - £ 2 (-Q)*W, k 2 ))~ {-Q) k ^U k , 10(^-1,^ + 1)) 

k=l 

N-l 

- E (-Q) k - 2 Mk 2 ,k)) (4.90) 

k=k 2 +2 

U \tt(h, N-l)) = (-Q)- N+2 U fea)) (4.91) 

where 2 < h < N - 3 and 3 < k 2 < N - 2. 

Following the same procedure of one-pseudoparticle case we again define new states 
in order to have consistency between bulk and boundaries terms. In addition to (|4.32|) 
we have the following new states 

U \Q(l,k 2 )) = \Q(0,k 2 )), U \n{k 1 ,N-l)) = \Q(ki,N)) 
U \Q(1,N-1)) = \n(0,N-l)) + \fl(l,N)) (4.92) 

Acting with H on these new states, we get 

H\n(o,k 2 )) = 2(Q + Q- 1 )\n(o,k 2 )) + \n(o,k 2 -i)) + \n(o,k 2 + i)) 

+ k 2 )) + (-Qf- 2 \n(k 2 , N - 1)) (4.93) 

H\n(h,N)) = 2(Q + Q- 1 )\Q(k 1 ,N)) + \Q(k 1 -l,N)) + \Q(k 1 + l,N)) 

+ \Q(k u N-l)) + (-Q)- N+2 \Q(1, h)) (4.94) 

Substituting these results into the eigenvalue equation, we get the following boundary- 
conditions 

A(k 2 , N + h) = (-Q) N - 2 A(k u k 2 ). (4.95) 
The parametrization for the wavefunctions 
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A(h, k 2 ) = A 12 && + Anfr& (4.96) 
solves the eigenvalues equations provided 

E 4s = 2(Q + Q- 1 ) + a + + 6 + (4-97) 

and the boundary conditions ( J4.95|) provided 

= (~Q) N - 2 ^ , & = (~Q) N - 2 ^ Z N = (-Q) 2{N - 2) (4.98) 

The phase shift produced by the interchange of two pseudo particle is again given by 

(El). 

We thus arrive to the Bethe ansatz equations which fix the values of £i and £2'- 



i+e+(Q+Q-% 



1 l + ^ + ^Q + Q-i)^/' 
= {-Qf {N - 2) (4-99) 

4.2.3 Two-pseudoparticles and impurities 

In the sectors 4s < r < 6s, we find states with two interacting particles and impurities. 
Let us now consider states with two pseudoparticles and one impurity, for instance, 
of type s — 1. We can recall the previous section to use all equations presented in 
the case of two-pseudoparticles and one impurity with periodic boundary conditions. 
Differences came only from boundary conditions terms 

Alfa, h, h + N) = (-Q) N - 2 Ai +1 (k u fc 2 , h) j = 1, 2, 3 (mod3) 

which read now 

(_Q)-iV+2£iV _ A 123 _ A 132 ,_Qy N+2 £N _ ^312 _ ^213 
^312 ^321 " ^231 ^132 

(-or™*? = 4^=4^, * N = = (-q) 2{n - 2) ( 4 - io °) 

^213 ^123 
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The eigenvalues of H are again given by ( |4.97| ) but with different Bethe equations, 
which have now a additional factor: 

" " [ ~ Q) JL\ i + + + ' 

ef- 4 = i = (-Q) 2M (4.ioi) 

It follows from our experience with the periodic boundary cases that these results 
can be extended to a generic sector. So, we will conclude this subsection with their 
generalization 

4.2.4 General eigenstates 

In the sector r = 2sp, we expect that the p-pseudoparticle phase shift will be a sum of 
p(p — l)/2 two-pseudoparticle phase shifts and the energy is given by 

e p(2s) = E {Q + Q- 1 + £n + C 1 } (4.102) 

n=l 

The corresponding eigenstates are 

*r(£i, 6, -4) = E fc 2,-, fep) k 2l (4.103) 

l</ci<...<fc p <A r -l 

where \Q(ki, k 2 , k p )) = ®f =1 |f2(/cj)) and the wavefunctions satisfy the following 
boundary conditions 

A(Aa, fc 2 ,.., k p , N + k 1 ) = (-Qf-^Aih, k 2 ..., k p ) (4.104) 

which produces an extra factor in the Bethe Ansatz equations. It is not all: in a sector 
r we may have pseudoparticles and impurities which play here the same role as in the 
periodic case. It means that for a sector r with I impurities with parameters £i, 
and p pseudoparticles with parameters &+ p the energy is given by ( |4.102| ), and 

the Bethe equations do not depend on impurity type and are given by 

with a = / + 1, / + 2, Z + p , p > 1, and 

e 2p te + i ■ ■ • = (-q) kjv ~ 2p+2) (4.106) 
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where f = £162 • • • • • • £ l+p . 

The additional factor in the Bethe Ansatz equations is a consequence of the nonlo- 
cality of Uq which generates new boundary conditions depending on the sector through 
the number of pseudoparticles and on the quantum supergroup parameter q via the 
relation Q + Q' 1 =Stry A (q 2p ). 

4.3 Free boundary conditions 

It is for free boundary conditions that the graded TL Hamiltonian naturally commutes 
with the quantum supergroup U q (X M \ n ). We expect that all procedures developed for 
the coordinate Bethe ansatz with free boundary conditions in |29) for the spin-1/2 
XXZ chain can be used here. To show this we recall the previous subsections, taking 
into account t.b. = 0, where almost all equations can be seized for the free boundary 
conditions eigenvalue problem. Indeed, we have to solve the bulk independently of the 
boundary terms. 

4.3.1 One-pseudoparticle eigenstates 

In this sector, the eigenstate is given by ( |4.79|) : 



where \Q{k)) is again given by ( |4.3| ). 

The action of H on the states \Q{k)) gives us the following eigenvalue equations 



<M0 = £ A(k) \n(k)) 



(4.107) 



k=l 



(E 2s -Q- Q^)A(k) — A(k - 1) +A(k+ 1), 2 < k < N - 2 



(4.108) 



At the boundaries, we get more two slightly different equations 



(E 2s - Q - Q- l )A(l) 
(E 28 -Q-Q- 1 )A(N -1) 



A{2) 

A(N -2) 



(4.109) 



We now try as a solution 



A(k) = A(0)£* - A(-0)£ 



-k 



(4.110) 
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where £ = e ld , 9 being the momenta. Substituting this in equation ( [4.1 08| ) we obtain the 
energy eigenvalue associated with a free pseudoparticle with free boundary conditions 

E 2s = Q + Q~ 1 +£ + r 1 (4-111) 



We want equations ( [4.108 ) to be valid for k — 1 and k = N — 1 also, where A(0) and 



A(N) are defined by (|4.110|) . Matching ( |4.108|) and ( [4.109| ) we get the end conditions 



A(0) = and A(N) = (4.112) 

implying that A(6) =A(—6) and £ 2Ar = 1, respectively. A(6) it now determined ( up 
to a factor that is invariant under 9 < — > —9), to be equal to £ _Af . 

4.3.2 One pseudoparticle and impurities 

Differently from the previous cases, due to the lack of periodicity, the impurity positions 
are fixed. So, they have a different role in the eigenvalue problem with free boundary 
conditions. For instance, let us consider the case of one impurity of the type s — 1, 
with parameter £i and one pseudoparticle with parameter £ 2 . For this case we obtain 
the following eigenvalue equations 

(E 2s+1 -Q-Q- 1 )A 1 (k 1 ,k 2 ) = A 1 (k 1 -l,k 2 )+A 1 (k 1 + l,k 2 ) 
{E 2s+x -Q-Q- l )A 2 {k u k 2 ) = A 2 (k u k 2 -1) + A 2 (k u k 2 + 1) (4.113) 

We also have two meeting conditions that arise because pseudoparticle and impurity 
may be neighbors (see (|4.20|) ) 

Ai(k, k) = A 2 (k } k + 2), A 2 (k + 1, k + 2) = A x (k, k + 1) (4.114) 

in addition to the two conditions to be satisfied at the free ends 

A 1 (k 1 ,N) = Q , A 2 (0,k 2 ) = (4.115) 

Now we try the following ansatz for the wavefunctions 

A 1 (k 1 ,k 2 ) = A 1 (e l ,e 2 )&&-A 1 (e 1 ,-e 2 )&& ha 

A 2 (h,k 2 ) = A 2 (e 1 ^ 2 )^^-A 2 (e 1 ,-0 2 )£ 2 - fci e i fe2 (4.116) 
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From ( |4.113| ) we get the energy eigenvalue 

E 2 s + i = Q + Q- 1 +& + & 1 (4.H7) 



and from ( 4.114 ) and (|4.115|) the following relations between the coefficients Ai 



A 1 {6 1 ,6 2 )& = A 1 (9 1 ,-9 2 )^ N , A 2 (9 1 ,9 2 )=A 2 (9 1 ,-9 2 ) 

A 1 (6 1 ,9 2 ) = A 2 (0i,0 a )£i , A 1 (9 1 ,-9 2 ) = A 2 (9 1 ,-9 2 )e i (4.118) 

from this we get 

e 2 N = i (4.H9) 

as the Bethe equation of ( (4.117] ). The coefficients Ai are determined up to a factor 
that is invariant under 9 2 < > —9 2 as: 

Ai(0i, 9 2 ) = e^ 2 N and A 2 (9 U 9 2 ) = & N . (4.120) 

In general, for the eigenstate with I impurities with parameters £i, and one pseu- 
doparticle with parameter which lies in a sector r, we can write 

*r(£i, =E E Mh, h+i)) \ (4.121) 

j=l I l<fci<...<fc !+1 <JV-l J 



The corresponding eigenvalue is given by Q4.111 ) , with £ = and the ansatz for 



the coefficients of the wavefunction becomes 



A j (9 1 ,...,e l+1 )= flRi Km (4-122) 



i=l 



Here we notice that the index j in the wavefunctions Aj(ki, ki+i) means that the 
pseudoparticle is at the position ki+ 2 -j- 

4.3.3 Two-pseudoparticle eigenstates 

For the sector r = 4s, beside eigenstates with impurities, we have an eigenstate with 
two pseudoparticles. We obtain the following eigenvalue equations 

{E As -2Q-2Q- 1 )A(k 1 ,k 2 ) = A{k 1 -l,k 2 )+A{k 1 + l,k 2 ) 

+A(ki, k 2 — l)+ A(ki, k 2 + 1) (4.123) 
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We have again two conditions to be satisfied at the ends of the chain 



A(0,k 2 ) = and A(k 1 ,N) = (4.124) 
In addition to this we have a meeting condition 

A(k, k) +A(k + l,k + l) + (Q + Q~ l )A(k, k + 1) = (4.125) 

Now we try the ansatz 

A(k u k 2 ) = A(9 1 ,9 2 )fr&-A(9 2 ,9 1 )&& 

- A{-0 1 ,6 a )£'*&+A{-0 a ,0 l )£ , »& , 419fi , 

- A(e 1 ,-e 2 )e l 1 ^ k2 + A(e 2 ,-e 1 )e l 2 G kl K ' J 
+ A(-e l7 -e 2 )^ 2 k > - A(-e 2 , -flOCr* 2 ^ 1 

Here we observe the permutations and negations of 6i and 6 2 . Substituting this ansatz 
m (14121 we obtain the energy eigenvalue for the sector with two pseudoparticles 



E 4s = 2Q + 2Q- 1 + & + + 6 + ^ (4-127) 
The ansatz ( |4.126| ) satisfy equations ( [4.124Q provided 



A(9 1 ,e 2 ) = A{-e u e 2 ) , A{6 2 ,e 1 ) = A{-e 2 ,e 1 ) 

A(9i, —9 2 ) = A{-9 u -9 2 ), A(9 2 ,-9 1 )=A(-9 2 ,-9 1 ) (4.128) 



and 



2N _ A(9 u -9 2 ) _ A(-9 1 ,-9 2 ) 2N _ A(9 2 ,-9 1 ) _ A(-9 2 ,-9 1 ) 
?2 " A(9 1 ,9 2 ) ^ A(-9 U 9 2 ) ' ^ " A(9 2 ,9 1 ) A(-0 2 A) 1 " ' 

Moreover, the meeting conditions are satisfied provided 

A(-0 1? -0 a ) A(0 a , 00 1 + Ci6 + (Q + Q _1 )6 



A(-0 2 ,-0O " A(^A) i + ^ + W + Q- 1 )^ 

a(9 1 ,-9 2 ) A(0 2 ,-eo _ i+er^ + CQ + Q- 1 )^ 



a(-0 2 ,0o " a(-^^ 2 ) i+er^ + ^+Q- 1 )^ 1 

Matching these conditions we get 



(4.130) 



27V _ B(-6 1 ,9 2 ) 2N _ B(-9 2 ,6i) 

4l ~ ' 42 " B(9 2 ,9 1 ) [AA6L) 

and 

A(0i, 2 ) = ^^(-^i, fc)^ 1 - (4-132) 
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Here we have used the usual free boundary notations 



PI 



B(9 a , 9 b ) 



9 h ) s(9 h , -9 n 



(4.133) 



where 



s(9 a ,9 b ) = i + Ub + (Q + Q- 1 )&- 



(4.134) 



Now let us consider the eigenstates with two pseudoparticle and impurities. The 
energy eigenvalue is the same of the two pseudoparticles pure state. The parame- 
ters associated with impurities are embraced in the definition of the coefficients of 
the wavefunctions. For instance, when we have an eigenstate of two pseudoparticles 
with parameters £2 and £3 and one impurity of parameter £1, the energy is given by 
fl£I27D and the Bethe equations by (|4.131|) , with & -> £ 3 and 9 X -»• 9 3 . But now the 
wavefunctions are different 



Ai (9 U 9 2 ,9 3 ) = {g}t 2 N & N B(-9 2 ,9 3 )& 
k 2 {9 u 9 2 ,9 3 ) = Bi-hMS 1 
A 3 (9 1 ,9 2 ,9 3 ) 



(4.135) 



where B(—9 2 ,9 3 ) is given by (|4.133| ) 
4.3.4 General eigenstates 

The generalization follows as in the previous cases. In a sector r with 
p pseudoparticles, we get 

p 



^r = E[Q+g- 1 +en+c 1 



n=l 



and the Bethe equations 



2N 



l+p 

n 



B(—9 a , 9b) 



b^a=l+l B(9 a ,9 b ) 



a = l,2,...,p 



The corresponding eigenfunction can be written as 



(4.136) 



(4.137) 



*r(fi, -,Q = J2 A(k u k 2 ,...,k l+P )\tt(k u k 2 ,...,k p )) 

fci<— <k t+p 



(4.138) 
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with 

A(h, k 2 , M = E £ p MOi, 82, Op) ■ ■ ■ & (4-139) 

p 

where the sum extends over all permutations and negations of 9i, 9 P and Bp changes 
sign at each such interchange. The coefficients in the wavefunction are given by 

A(9 h 9 2 , 9 P ) = f[Zj N II B (- e 3> Oi)^ 1 ( 4 - 14 0) 
j=i i+i<j<i<i+ P 



where B(—9j,9i) are defined in ( |4.133| ). 

For a sector r with I impurities with parameters £1, and p pseudoparticles 
with parameters £z + i, ...,£j +p the energy is given by ( |4.136| ) and the Bethe equations 



by ( [4.137 ) . Only the coefficients of the wave functions are modified 

Aj(0i, 02, 9 l+p ) = ■ ■ ■ 6)A(^+i, 2 , 9 l+p ). (4.141) 

The functions A, (£i£ 2 • • • 6) = ^i 1 ^ 2 ' ' ' where the index j characterizes the possible 
configurations of I impurities relative to the p pseudoparticles. Here a« are numbers 
which depend on the position of corresponding impurity relative to the pseudoparticles. 

We notice again that these results are valid for all TL Hamiltonians defined as 
projector of spin zero on the representations of the quantum supergroups U q (XM\ n ), 
characterized by the values of Q + Q" 1 = Try A (g 2p ). 

5 Conclusion 

Via the coordinate Bethe Ansatz, we obtained the spectra of a series of "spin " Hamil- 
tonians associated with the representations of the graded TL algebra arising from 
orthosymplectic quantum supergroups. We believe that the Bethe Ansatz technique 
used in this paper can be applied to all Temperley-Lieb models. 

Due to the TL operator being a one- dimensional projector, there is a linear com- 
bination of eigenstates of S z = diag(s, —s) , which is proportional to the singlet. 
In terms of this Bethe vector we get a unified treatment for all Hamiltonians, which 
shows the TL equivalence at Bethe Ansatz level. We find that all models have equiv- 
alent spectra, i.e., they differ at most in their degeneracies. Moreover, for the closed 
boundary conditions, the spectra of the lower- dimensional representations are entirely 
contained in the higher-dimensional ones. 
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Here we notice that this spectrum equivalence is, of course, a consequence of the 
TL algebra. Nevertheless there is in the literature a large class of Hamiltonians which 
are not derived from representations of the TL algebra which share the same property. 
The authors of reference [3(| developed a technique for construction of spin chain 
Hamiltonians which affine quantum group symmetry whose spectra coincides with the 
spectra of spin chain Hamiltonians which have non-affine quantum group symmetry. 

There are several issues left for future work. In particular, one can derive the 
partition functions in the finite size scaling limit to find the operator content for the 
systems constructed from these quantum chains. 
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